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We have investigated the classical phase dynamics of underdamped ferromagnetic Josephson junc- 
tions by measuring the switching probability both in the stationary and non-stationary regime down 
to 350 mK. We found the escape temperature to be the bath temperature, with no evidence of ad- 
ditional spin noise. In the non-stationary regime, we have performed a pump-probe experiment on 
the Josephson phase by increasing the frequency of the junction current bias. We show that an 
incomplete energy relaxation leads to dynamical phase bifurcation. Bifurcation manifests itself as 
premature switching, resulting in a bimodal switching distribution. We directly measure the phase 
relaxation time t v by following the evolution of the bimodal switching distribution when varying 
the bias frequency. Numerical simulations account for the experimental values of t v . 

PACS numbers: 74.50.+r,74.45.+j, 85.25.Cp 



The maximum dissipationless current that a Joseph- 
son junction can carry, J s , is not only given by the cou- 
pling energy between the two superconductors forming 
the junction, but also by the dynamics of their phase 
difference tp [l[. Thermal and quantum fluctuations of 
the phase, for instance, lower I s from the critical cur- 
rent value Ic = 7rA/2ei? n , where A is the superconduct- 
ing gap and R n the junction resistance in the normal 
state [2]. Furthermore, they introduce in a current-biased 
junction a probability distribution of the threshold cur- 
rent to switch into the dissipative state Two time 
scales are pertinent: the inverse plasma frequency u)q , 
and the phase relaxation time, r vi which is defined by 
the junction's damping Since the plasma frequency 
is in the gigahertz range, it is actually the microwave 
impedance, including that of the junction environment, 
that sets the relaxation. However, in large-area tunnel 
junctions, the effect of the environment is negligible as 
the quasiparticle resistance is smaller than the vacuum 
impedance and the phase dynamics is "intrinsic" in the 
sense that it only depends on the junction parameters 
1- 

In this Letter, we address the classical phase dynam- 
ics of large-area underdamped ferromagnetic Josephson 
junctions. We show that in the stationary regime, i.e. 
when the current bias frequency Ub is much smaller than 
the inverse phase relaxation time t" 1 , the mean switch- 
ing current (I s ) and its standard deviation a are as ex- 
pected for thermal fluctuations. Whereas in the non- 
stationary limit, i.e. when the bias frequency is com- 
parable to the inverse phase relaxation time, (7 S ) has a 
finite probability to jump to zero, thus showing an ap- 
parent zero critical current. The premature switching 
reveals the phase bifurcation. The exponential rise of 
the premature switching probability N\ as a function of 
the bias frequency u>b gives the phase relaxation time t v 
in agreement with the numerical simulations. Moreover, 
no effects of the magnetization on the phase dynamics at 



zero voltage are observed. 

The phase dynamics of a Josephson junction is usu- 
ally described within the Resistively Shunted Josephson 
junction (RSJ) model [B| as 

tp + f3(p + simp = rj b (t), (1) 

where (3 = (RCujq)^ 1 , luq = d^elcJfrC is the plasma 
frequency and rjb = Ib/Ic, where is the time de- 

pendent bias current. In our experiment it is a sawtooth 
of frequency wt>. The time is normalized to r = u>ot, and 
the dot denotes d/dr. Here C is the junction capacitance 
and R the quasiparticle tunneling resistance. Damping /3 
is set by the RC product. Escape by thermal fluctuations 
has been largely investigated since the pioneering work 
of Fulton and Dunkelberger Q . Thermal escape is usu- 
ally described using the Kramers theory in the harmonic 
approximation [?J ■ The kinetic energy of the order of the 
Josephson energy, and much higher than the thermal en- 
ergy, makes the harmonic approximation incorrect. The 
anharmonicity produces bifurcation Q in the phase dy- 
namics that in the end strongly modifies the switching 
distribution probability, rendering it bimodal. 

Bifurcation usually occurs when an anharmonic oscil- 
lator is driven above a certain critical amplitude close to 
the resonant frequency. This creates two stable attractors 
in the phase space and hence two macroscopically differ- 
ent states. An example of a system in which bifurcation 
occurs is the Duffing oscillator Q. Dynamically driven 
bifurcation is often referred to as a kinetic phase transi- 
tion The noise allows for stochastic transitions be- 
tween the two attractors leading to the finite probability 
of finding the system in each state. Microwave-induced 
bifurcation for the Josephson junctions has been observed 
in the stationary regime [ll[ and has been realized exper- 
imentally in the pioneering work by Yurke et al. [l2| for 
parametric amplification 3, 3|- 

In the non-stationary regime [HI], it is the initial ki- 
netic energy that induces the phase bifurcation. Two 
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phase states, shown in Fig. 1(a), are achieved: Either 
the phase relaxes to the bottom of the potential well re- 
sulting in a switching current / s equal to that in equilib- 
rium, or the phase does not relax and the escape occurs 
for much smaller current bias, of the order of the retrap- 
ping current I r . In strongly underdamped junctions as 
these, the retrapping current is vanishingly small. The 
two solutions are illustrated in the Fig. 1(b), showing the 
numerically obtained V{t) from the RSJ model [Eq. (JTJ], 
considering a periodic bias current with a sawtooth wave- 
form. After the first escape at I s , the kinetic energy is 
large enough for premature escape at I r when the cur- 
rent is ramped up again. Thus in the non-stationary limit 
the kinetic energy induces a phase transition between the 
locked and the running state [lfjj]. This is illustrated in 
the calculated phase diagram (tp,(p) [Fig. 1(c)]. The full 
and dotted curves correspond to the locked and running 
state respectively. The initial kinetic energy, i.e., initial 
speed ip{t = 0) at the start of each cycle, determines the 
phase trajectory, and the dashed line is (approximately) 
the separatrix between the running and the locked state 
trajectories. 
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FIG. 1: (color online), (a) The switching mechanism in the 
non-stationary regime, where the phase can either relax to 
the bottom of the potential well and then escape (solid cir- 
cle), or stay near the top and escape early (open circle), (b) 
The calculated voltage V(t) across the junction (full curve), 
with the imposed current sawtooth ramp (dotted curve). Ei- 
ther the switching takes place as in equilibrium, around I s , 
or around the retrapping current I r . (c) The phase diagram, 
calculated for f3 — O.Olwo, uJb = O.OOIojq and r/b = 1.2. The 
attractor (solid curve) corresponds to the locked state, while 
the dotted curves correspond to the running state. In the case 
of the early escape, the phase trajectory is near the separatrix 
(dashed curve), (d) Numerically obtained distribution of the 
initial phase speed at the start of each cycle (ip(0)), as a func- 
tion of the ramp frequency. The error bars are augmented 
5 times. Full line is a linear fit, and the dashed line is the 
bifurcation threshold corresponding to the separatrix. 

The dynamical bifurcation depends on the bias fre- 



quency, bias amplitude and damping. It turns out that 
it occurs only when the phase relaxation time, r v , is com- 
parable with the driving force period 2ir/ojb- Note that 
t v in strongly underdamped junctions is 4 orders of mag- 
nitude larger than the resonant period 2tt /luq- Therefore, 
the out-of-equilibrium conditions shift the time scale for 
bifurcation from the resonant frequency, as expected in 
the stationary regime, down to the phase relaxation time. 
The numerical solution of the RSJ model shows that no 
premature switching is expected when t v <§C 2ir/u>b, while 
it is indeed observed at t v ~ 2ir /ujb- For t v 3> 2tt /u>b we 
find only the premature switching solution. Physically, in 
this case the phase has not relaxed when a small driving 
force (i.e., a current bias equal to the retrapping current) 
allows it to escape. The effect of noise in the numerical 
model is replaced by the fluctuations of the initial ki- 
netic energy for each current sweep due to an incomplete 
relaxation. These fluctuations allow for stochastic tran- 
sitions between the locked and the running state. The 
numerically obtained initial speed (averaged over many 
periods) {(fi(0)) as a function of the ramp frequency in 
Fig. 1(d) sums up the effect of the bias frequency on the 
kinetic phase transition. The fluctuations of the initial 
speed are given by the error bars. The dashed line is the 
separatrix, so when the initial speed is over the bar, there 
is a finite probability for bifurcation. 

In the nonstationary limit, the probability for prema- 
ture switching TVi , obtained by numerically solving the 
RSJ model 17] [Eq. (fTJ] for a sawtooth ramp of ampli- 
tude r\b and frequency ojb, is well described by 



Ni = 1 — Aexp(— T v ujb) = 1 — 1.8 exp 



0.76 



Vb^b 



(2) 



At t v ~ 2n/uJb 1 the current switching distribution is 
bimodal, and Ni is the surface of the histogram at 
I r . In a pump-probe experiment, the relaxation time 
is usually obtained by tracking the time needed to reach 
equilibrium after an external perturbation is turned off. 
Here, similarly, the phase relaxation time can be directly 
probed by measuring the early switching probability N% 
as a function of the bias frequency. This offers an alter- 
native way to measure oscillation damping [l8| . 

Experimental. We fabricate large-area superconduc- 
tor-insulator-ferromagnet-superconductor (SIFS) tunnel 
junctions by electron gun evaporation through shadow 
masks in ultrahigh vacuum. The superconductor is Nb, 
the insulator is AI2O3, the ferromagnet a PdNi alloy with 
10 % Ni, and junction surface is 0.6 x 0.8 mm 2 . The fabri- 
cation technique is described elsewhere [lj|. In this Let- 
ter, we investigate escape in Josephson junctions showing 
only 7r coupling and corresponding to a thickness of the 
PdNi ferromagnetic layer dp varying from 70 to 100 A. 
We concentrate on the n junction as it has been used as a 
phase battery in superconducting circuits 2^, 21 1 , and it 
has been argued that it may be used in superconducting 
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Q-bits [22| . The coupling between the phase and the spin 
dynamics of the PdNi thin layer is an important issue. 
Furthermore, from a technical point of view, the ferro- 
magnetic layer suppressed the critical current by a factor 
of 100, reducing substantially the power dissipated upon 
switching, keeping the heating negligible. 

Junctions have been measured in a 3 He cryostat from 
4.2 K down to 350 mK. The current-voltage curve of a 
typical SIFS junction at 350 mK is presented in Fig. 2(a). 
The curve is that expected for a strongly underdamped 
junction. The critical current is Ic — 130 \xA. The junc- 
tion resistance in the normal state is R n = 0.2 ft, and 
capacitance is C = 42 nF, as estimated from the AI2O3 
dielectric constant and the Fiske steps. We obtain the 
plasma frequency u>q = 3 x 10 9 rad/s, and the quality fac- 
tor Q = = 3000. We estimate the retrapping current 
to be I r ~ AIc/ttQ ~ 50 nA. Finally, the power dissipated 
at switching is of the order of 1 nW. As the junction 
volume is about 10 5 /im 3 , the cooling power from the 
electron-phonon coupling leads to an irrelevant increase 
of the electron bath temperature. 

The switching current distribution is measured using a 
counter. The junction is current-biased with a sawtooth 
waveform, and the switching current is easily obtained 
by multiplying the time before dissipation and the speed 
of the ramp. Positive and negative bias give the same 
result. All of the incoming and outgoing lines are filtered 
with band pass filters with the cutoff frequency of 1 MHz, 
at 4.2 K and room temperature. 
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FIG. 2: (color online), (a) Typical IV curve at 350 mK. 
(b) Distribution width as a function of temperature for junc- 
tions with 70, 80, and 85 A of PdNi (triangles, circles, and 
squares, respectively) . The solid curve is a fit for the junction 
with 80 A in the low damping regime, (c) The measured his- 
tograms P(J) (dotted curves) for T = 0.5, 0.8, 1.1, 2.0, 3.0, 
and 4.2 K, respectively, from the right, and a fit (solid curves), 
for a junction with 80 A of PdNi. The effective temperature, 
obtained as a fitting parameter, is almost equal to the bath 
temperature. 



Thermal Escape. The switching current distribution 
of the junction with 80 A of PdNi, measured at differ- 
ent temperatures, is presented in Fig. 2(c) for the bias 
frequency of 37 Hz. This time scale is much longer than 
the phase relaxation time, and hence the phase follows 
adiabatically the drive. The mean switching current in- 
creases linearly with decreasing temperature jljjj . The 
width narrows as expected from escape induced by ther- 
mal fluctuations. The fit (solid curves) is obtained using 
the Kramers theory 0, [23| for low damping. 

The distribution width a as a function of bath tem- 
perature T is presented in Fig. 2(b) for three different 
junctions, with d F = 70, 80 and 85 A of PdNi. The 
solid curve is the fit in the low damping regime, which 
takes into account the temperature-dependent quasipar- 
ticle resistance R No phase noise induced by spin 

noise is observed. Note that here we focus on the clas- 
sical phase dynamics, and no crossover to the quantum 
tunneling is expected down to 10 mK. Spin noise orig- 
inating from thermal fluctuations of the magnetization 
can affect the escape through the fluctuations of the crit- 
ical current SIm or through the magnetization-induced 
phase noise 5<pm- The critical current depends on the ex- 
change energy as 7c~e x p( — <£f/£f) cos(cZf/£f) |25] , with 
£f = y hD / E ex , where E ex is the exchange energy and 
D the diffusion coefficient. For D ~ 10cm 2 /s, E ex ~ 25 
mV and dp ~ 80 A, we have c?f/£f ~ and SIm/Ic 
goes as (SE ex /E ex ) 2 . Therefore, for SE ex ~ fc^T, we get 
SI m ~ 10 nA, which is below our resolution. The over- 
all SifM produced by thermal domain wall motion and 
intradomain oscillations of the magnetization is small in 
large-area junctions. Thermal phase fluctuations account 
well for the switching current distribution measured at 
equilibrium. 

Phase Relaxation Time. In order to investigate the 
change in the switching current distribution when the 
phase is nonstationary, we systematically increase the 
bias frequency. We observe that current switching distri- 
bution changes drastically and becomes bimodal above 
5 kHz. With increasing bias frequency, the number of 
counts around I s decreases, while the number of counts 
around I r increases. To increase sensitivity, instead of 
the retrapping current I r , we measured the current cor- 
responding to the first Fiske resonance If ~ 15 /iA (2(|. 
This is shown in Fig. 3(a), where we report the switching 
current distribution for three different bias frequencies. 
The total number of sweeps is 2000. Above 20 kHz, the 
mode around I s disappears. The histograms at high bias 
frequency show only a tiny reduction in the mean value 
of / s with respect to that measured at equilibrium, while 
almost no change in their width indicates negligible ther- 
mal effects. It is worthwhile to mention that we have ob- 
served this behavior in Josephson junctions without any 
ferromagnetic layer when the critical current is almost 
canceled by applying a quantum of magnetic flux in the 
junction (i.e., around the first minima in the Fraunhofer 



4 



140 




FIG. 3: (color online), (a) Switching histogram of the junc- 
tion with 85 A of PdNi taken at 350 mK with the sweep 
frequency of 4 kHz (solid curve), 6 kHz (dotted curve) and 
12 kHz (dashed-dotted curve), (b) Early switching probabil- 
ity Ni as a function of ramp frequency, measured at 350 mK 
for the junctions with 85 (open circles), 100 (solid squares) 
and 90 A (crossed squares) of PdNi. The dotted lines are a 
theoretical prediction [Eq. (2)]. The solid curves are obtained 
from the same equation, when keeping t v as the fitting pa- 
rameter. Inset: The dependence of the bifurcation onset fre- 
quency v* of temperature, measured on a junction with 85 A. 
The solid line is a fit taking into account the measured quasi- 
particle resistance R. 



pattern). 

The premature switching probability N\ as a function 
of bias frequency is presented in Fig. 3(b). The dot- 
ted curves are the solutions of the RSJ model [Eq. (2)] 
without fitting parameters, while the solid curves are the 
best fit, giving t v = 157, 207, and 282 fia, respectively, 
for junctions with 85, 100, and 90 A of PdNi [27j]. The 
discrepancy could arise from the fact that in the numer- 
ical model we do not account for thermal noise. The 
frequency at which the bifurcation starts, v* , is repro- 
duced numerically. An independent way to measure the 
RC product is to measure the delay of V(t) relative to 
the start of the ramp, and the obtained values are in 
accord with the experimental t v . 

Finally, we show the measured bias frequency at which 
bifurcation starts, v* , as a function of temperature [in- 
set in Fig. 3(b)]. The fit is obtained numerically as a 
solution of Ni(v*) — Q, taking into account the measured 
quasiparticle resistance. 

Conclusions. We measure the switching histograms 
of the strongly underdamped SIFS Josephson junctions. 
In the stationary regime, the escape temperature is the 
bath temperature indicating negligible magnetic-induced 



noise. Therefore at least the large-area 7r junction seems 
to provide a noiseless phase battery. In the nonstation- 
ary regime, when the junction is ramped before complete 
phase relaxation, we observed a kinetic phase transition 
induced by phase bifurcation. The probability of measur- 
ing a finite switching current decays as a function of the 
bias frequency, allowing to measure the phase relaxation 
time as in a pump-probe experiment. From a practical 
point of view, dissipationless current of a Josephson junc- 
tion can be the turned on and off by changing the bias 
frequency. 

We thank J. Gabelli, B. Reulet, T. Kontos and Z. 
Radovic for stimulating discussions. 
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